Vacancy supersolid of hard-core bosons on the square lattice 
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The ground state of hard-core bosons on the square lattice with nearest and next-nearest neigh- 
bor repulsion is studied by Quantum Monte Carlo simulations. A supersolid phase with vacancy 
condensation and 'star' diagonal ordering is found for filling p < 0.25. At fillings p > 0.25 a su- 
persolid phase exists between the star and the stripe crystal at p=0.5. No supersolid phase occurs 
for p > 0.25 if the ground state at half-filling is either a checkerboard crystal or a superfluid. No 
commensurate supersolid phase is observed. 
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Experimental advances in the field of ultracold atoms 
in optical lattices [l[ have given renewed impetus to the 
investigation of novel phases of matters, especially those 
displaying simultaneously different types of order. One 
such phase is the supersolid, featuring both diagonal and 
off-diagonal long range order, which has been the subject 
of much theoretical speculation Q ■ Recent claims of ob- 
servation of this phase in solid helium 0] have generated 
some controversy 

Theoretical studies have yielded strong evidence of su- 
persolid phases of lattice bosons, for various types of 
model interactions among atoms, as well as of lattice ge- 
ometries [H @, 0, 13 ■ In the presence of nearest- neighbor 
particle hopping and on-site hard core repulsion, the su- 
persolid phase is only observed on the interstitial side of a 
commensurate solidphase (e.g., for 1/3 < p < 2/3 on the 
triangular lattice) [8|. In contrast, doping with vacancies 
results in the coexistence of an insulating crystal and a 
superfluid by the formation of a domain wall Q . Phase 
separation of vacancies is also observed in ah initio sim- 
ulations of helium crystals [Tcj. This seems remarkable, 
as Bose condensation of vacancies has long been regarded 
as the paradigm for supersolidity Q ■ The purpose of this 
work is to gain understanding in the asymmetry between 
the behavior of vacancies and interstitials, and explore 
physical conditions that underlie a vacancy supersolid 
phase in lattice bosons. 

Our starting point is the well-known Hamiltonian 

h = — tj2(Hh + h - c -) + y i hifi i 

(y) (y) 

A square lattice with periodic boundary conditions of 
N = L x L sites is assumed. The boson density (filling) 
is p = Ng/N, where the number of particles Nb is deter- 
mined by the chemical potential /i. The operator Sj(Sj) 
creates (annihilates) a hard-core boson on site i, with a 
maximum occupation number hi = feJSj of one particle 
per site. The first term of JT]) describes particle hopping 



to a nearest-neighboring site with amplitude t, which is 
our energy scale, t = 1. The second and third terms rep- 
resent repulsive interactions between bosons on nearest 
and next-nearest neighboring sites, respectively. 

Previous studies of the ground state of (p} have yielded 
evidence of three possible phases at half-filling ( p=0.5): 
a superfluid, a checkerboard solid and a stripe solid. The 
latter two are commensurate and insulating phases. Dop- 
ing the stripe crystal away from half filling yields a super- 
solid phase, whereas a first-order quantum phase tran- 
sition separates the checkerboard crystal from a super- 
fluid For sufficiently strong next-nearest neighbor 
repulsion, the Hamiltonian (|TJ) also features a commen- 
surate crystal phase at quarter filling, referred to as 
"star" (Fig. [T]). It has also been shown that a super- 
solid phase can occur, upon doping the star solid above 
quarter-filling Q. 
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(a) (b) 

Figure 1: (Color online). Classical, degenerate star crystal 
phases at quarter-filling on the square lattice. Filled circles 
represent particles, empty circles empty sites. Quantum fluc- 
tuations lift the degeneracy, selecting (a) as the ground state 
if Vi < 2V2, (6) otherwise. 

In this work, we systematically investigated the ground 
state phase diagram of ([TJ) . Our main findings are: 

i) Doping a star solid with vacancies always gives rise to 
a supersolid phase. 

ii) Between a star and a stripe solid a supersolid phase 
can be observed; phase transitions are continuous. 

Hi) If the phase at half filling is not a striped solid but 
a checkerboard solid or a superfluid, no supersolid phase 
exists. First order phase transitions separate superfluid 
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we find no commensurate 



from crystal phase (s). 
iv) In agreement with Rcf. 
supersolid phase [l2T |. 

Our Monte Carlo simulations are based on the Worm 
Algorithm in the lattice path- integral representation [l3| , 
using the implementation described in Ref. Q- The sim- 
ulations are at low temperature (typically j3 = L), in or- 
der to extract the properties of the system in its ground 
state. 
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Figure 2: (Color online). Ground state phase diagram of {T} . 
Lines are guides to the eye. Symbols lie at computed phase 
boundaries, namely between a) a "star" crystal and a super- 
fluid at p = 0.25 (open circles) b) a "striped" (n, 0), (0, tt) crys- 
tal and a superfluid at p — 0.5 (filled boxes) and c) a super- 
fluid a "checkerboard" crystal (filled circles), also at p = 0.5. 
Statistical errors are smaller than symbol sizes. Results at 
p = 0.5 are consistent with previous works [jj. 



In order to characterize the various phases, we com- 
puted the superfluid fraction p$ using the well-known 
winding number estimator 15j , as well as the static struc- 
ture factor 
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Figure 3: (Color online). Density p versus chemical potential 
p for Vi = 8, V% = 3.5. Statistical errors are smaller than 
symbol sizes. A first-order quantum phase transition sepa- 
rates the "star" crystal at p = 0.25 and the "checkerboard" 
crystal at p = 0.5 from a superfluid phase at intermediate 
densities. Results shown are for a square lattice with L=24. 
A continuous phase transition to a vacancy supersolid occurs 
below p = 0.25. 




S(Q) = 




where (. . .) stands for thermal average. Presence of 
crystalline long-range order is signaled by a finite value 
of i5(Q) for some specific wave vector. In particular, 
Q = (tt, tt) is the wave vector associated to checkerboard 
order at half filling, whereas Q = (n, 0), (0, ir) signals 
striped order at half filling, as well as star order at quar- 
ter filling. 

Fig. [2] summarizes our results for the ground state 
phase diagram of (JTJ) . Three distinct phase boundaries 
are identified. The first one (open circles) separates a su- 
perfluid from a star crystal at quarter filling. In this 
work, we focus our attention on the region above the 
open circles, i.e., wherein a star crystal exists. The other 
two phase boundaries refer to the behavior of the system 
at half-filling, separating a stripe crystal (filled boxes) 



Figure 4: (Color online). Superfluid density ps (upper panel) 
and static structure factor 5(Q) (lower panel) for Vi = 8, 
V2 = 3.5 for two system sizes and the two wave vectors Q = 
(tt,tt) and Q = (n,0), (0,n). Statistical errors are smaller 
than symbol sizes. 

[checkerboard crystal (filled circles)] from a superfluid. 
We will now discuss the physical behavior at other den- 
sities as a function on V± and V%. 

Fig. [3] shows the density computed as a function of 
the chemical potential for V\ = 8, V2 — 3.5. The ground 
state of the system is a checkerboard crystal at half fill- 
ing. Density jumps signal first-order phase transitions 
above quarter- and below half-filling; everywhere else the 
curve is continuous albeit with a clear discontinuity of the 
first derivative (i.e., of the compressibility) on approach- 
ing p = 0.25 from below. Fig. 2] displays the superfluid 
density as well as the static structure factor, computed 
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Figure 5: (Color online). Density p versus chemical potential 
p for Vi — 6.5, V2 = 3.5 and linear system size L = 24. 
Statistical errors are smaller than symbol sizes. A continuous 
quantum phase transition separates the star crystal from a 
supersolid on the vacancy side ( p < 0.25), whereas a first- 
order phase transition from the stripe crystals to a superfiuid 
exists on the interstitial (p > 0.25) side. 

for both (7r,0) and (ir,n). The superfiuid density is ev- 
erywhere finite, except at commensurate fillings, whereas 
the static structure factor is finite at and near commen- 
suration. 




Figure 6: (Color online). Superfiuid density ps (upper panel) 
and static structure factor S(Q = (tt, 0), (0, tt)) (lower panel) 
for Vi = 6.5, V% — 3.5 for two system sizes. System size 
is L — 24. A supersolid phase exists only on the vacancy 
side below quarter-filling. The ground state of the system at 
half-filling is superfiuid. 

Based on these results, we conclude that the system 
is superfiuid for 0.25 < p < 0.5, with first-order quan- 
tum phase transitions to the star (checkerboard) crystal 
at p = 0.25 (p = 0.5). For < p < 0.25, on the other 
hand, the system features two continuous phase transi- 
tions, one from a superfiuid to a vacancy-based supersolid, 
corresponding to the change in slope of the p{p) curve, 
and the other from the supersolid to the star crystal at 




Figure 7: (Color online.) Density p versus chemical poten- 
tial p for Vi = 5, V2 = 5 and linear system size L = 24. 
Statistical errors are smaller than symbol sizes. Continuous 
quantum phase transitions separate both star and stripe crys- 
tals from supersolid phases on both sides. A continuous phase 
transition from a supersolid to a superfiuid takes place below 
quarter-filling. 

p = 0.25. It is interesting to compare the results right 
above and below quarter filling (Fig. [4j in order to appre- 
ciate the effectiveness in distinguishing a homogeneous 
supersolid phase from coexistence of two phases. 

A similar scenario is observed for (V\,V2) pairs for 
which the ground state of the system at half-filling is 
superfiuid, where the only difference is the absence of 
the crystalline phase at p = 0.5. In this case too, a 
vacancy-based supersolid exists below quarter-filling and 
a first-order phase transition separates the star crystal 
from a superfiuid above quarter-filling. The superfiuid 
phase extends all the way to half-filling as the competi- 
tion between nearest-neighbor and next-nearest-neighbor 
interactions does not result in any crystal ordering. Re- 
sults for one example in this region, i.e., V\ = 6.5, V2=3.5, 
are shown in Figs. [5] and [6] 

This behavior should be contrasted to the one shown 
in Figs. [7] and [8] for the case V\ = V2 = 5. This pa- 
rameter choice stabilizes a star crystal at quarter-filling 
and a stripe crystal at half-filling. A supersolid phase 
exists everywhere between the star and the stripe crys- 
tal. A vacancy-rich supersolid phase occurs below filling 
0.25, with a continuous phase transition to a superfiuid at 
lower densities. A phase transition also occurs in this case 
above quarter filling, between a star and a striped super- 
solid; its location is signaled by a "kink" in the numerical 
value of the static structure factor (Fig. [8j lower panel). 
This phase transition was claimed to be of first order in 
Ref. 1|| based on a discontinuous jump in the value of the 
quantity 1 5* (7r, 0) — S(0,ir)\ across the transition. If the 
values of V\ and V2 are such that a checkerboard crystal 
or a superfiuid emerges at half-filling no supersolid phase 
intervenes in the 0.25 < p < 0.5 interval. 

The physical mechanism underlying the presence of a 



4 




ranged (e.g., dipolar) interactions among particles. The 
ratio between V\ and V2 for polar molecules is 2-^(2). 
The star solid, and the corresponding vacancy supersolid 
mig ht thus be seen in experiment for Vi > 10 (sec Fig. 2) 
|17| . albeit with rather strong requirements on the den- 
sity and temperature. Observation of a supersolid seems 
thus possible with existing techniques. 

This work was supported in part by the Natural Sci- 
ence and Engineering Research Council of Canada under 
research grant G121210893, by the Alberta Informatics 
Circle of Research Excellence and by the Swiss National 
Science Foundation. Part of the simulations ran on the 
Hrcidar cluster at ETH Zurich. 



Figure 8: (Color online). Superfluid density ps (upper panel) 
and static structure factor S(Q = (71", 0), (0, 7r)) (lower panel) 
for Vi = 5, Vi = 5. System size is L — 24. A supersolid phase 
exists everywhere between quarter- and half-filling, and below 
quarter-filling. 

supersolid phase in this model away from commensurate 
fillings, even on lattices other than the square, is whether 
or not defects such as interstitials or vacancies can move 
without frustration. On the triangular lattice at filling 
p—1/3, and in the presence of nearest-neighbor repulsion 
only, interstitial particles can freely hop around, hence 
Bose-condensing. At a finite concentration of vacancies, 
on the other hand, the system can lower its energy by 
forming a domain wall, i.e., vacancies phase separate [8|. 
The roles are reversed in our case, on the square lattice 
at quarter-filling. There is no cost in moving a vacancy 
around, in either the Fig. [lja) or[lj6) scenarios. Vacan- 
cies will thus condense. Conversely, the lowest energy 
needed to create an interstitial in the case of Fig. QJ6) 
(i.e., V\ > 2V2) is 4V2, which requires that a crystal atom 
also move to an empty row. The energy can only be low- 
ered by higher order processes in t/Vi, which gives rise 
to the same domain wall argument invoked above for va- 
cancies. This is the phase separation shown in Fig. [3J 

The lowest energy needed to create an interstitial in 
the case shown in Fig. QJa) (i.e., for V\ < 2V2 ) is 2Vi, 
which corresponds to an interstitial between two crystal 
atoms. Those two crystal atoms can now hop around 
at no potential energy cost, which corresponds to the 
supersolidity seen in Fig. [7] This simple perturbative 
argument fails to account for the scenario of Fig. [5j as 
higher order processes renormalize the transition point. 
The same physical considerations also easily explain the 
lack of a supersolid phase of either vacancies or intersti- 
tials on the kagome or honeycomb lattices (l6j . 

In an optical lattice, interactions among nearest and 
next-nearest neighboring atoms could arise from long- 
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